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We study the factorization of the dijet cross section in e^e“ annihilation using the generalized 
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Cambridge/Aachen jet algorithms as special cases. In order to probe the characteristics of the jet 
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a < 2. The factorized jet function and the soft function are well defined and infrared safe for all 
the jet algorithms except the kT algorithm. The kT algorithm (a = 2) breaks the factorization 
since the jet and the soft functions are infrared divergent and are not defined for a = 2, though 
the dijet cross section is infrared finite. In the jet algorithms which enable factorization, we give a 
phenomenological analysis using the resummed and the fixed-order results. 
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I. INTRODUCTION 


In high-energy scattering, strong interaction plays a signihcant role both in the high- 
energy region and in the low-energy region. The scattering of energetic quarks and gluons 
with radiative corrections occur in the high-energy regime, in which perturbative QCD can 
describe the effects of the strong interaction. On the other hand, strong interaction is also 
responsible for pulling the quarks and gluons out of the initial hadrons through the parton 
distribution functions, and for the hadronization of the energetic scattered partons through 
the fragmentation functions. These are nonperturbative aspects of the strong interaction 
in the low-energy regime. Furthermore the strong interaction affects all these processes at 
different energy scales at the same time, therefore the effects of the strong interaction are 
entangled in a complicated way all through the scattering process from beginning to end. 

It is important to disentangle the short- and the long-distance effects of the strong inter¬ 
action to have theoretical predictive power. This is achieved by the factorization theorem in 
which the scattering cross section is factorized into the hard, collinear and soft parts. The 
factorization proof in various hard scattering processes has been a long-standing problem in 
QCD [HE]. Recently the advent of the soft-collinear effective theory (SCET) [3H5] facilitates 
the proof in a straightforward way. SCET is formulated such that the collinear and the soft 
modes are decoupled at the Lagrangian level, and the factorization of the scattering cross 
section into the hard, collinear and soft parts is intuitively transparent. 

In high-energy processes, hnal-state hadrons are typically clustered in collimated beams, 
which are called jets. Jets are important tools to exploit information on Standard Model or 
beyond. There are many jet dehnitions depending on the kinematics, or the detector design, 
etc. [6]. In this paper we will conhne ourselves to the study of the various jet dehnitions used 
in annihilation, but it can be extended to ep scattering or to pp scattering. The jet 
dehnitions are realized by providing the jet algorithms which give rules on how to combine 
particles in forming a jet. The jet algorithms include the Sterman-Weinberg algorithm [7], 
the cone-type algorithm [8], the JADE algorithm [9], the kx algorithm [101, the anti-fc^ 
algorithm m, the Cambridge/Aachen algorithm [T^ 113] to name a few. 

The basic issue in jet algorithms is how to combine adjacent hnal-state particles into a jet 
without ambiguity. The important ingredients of appropriate jet algorithms are such that 
they should yield infrared safety, and they can be employed conveniently both in theory and 
experiment. However, the implementation of a convenient jet algorithm in experiment is 
not always convenient in theory, or vice versa. 

The factorized form of the dijet cross section in e+e“ annihilation can be obtained in the 
framework of SCET. The cross section is factorized into the hard, collinear and soft parts. 
However, the factorization does not hold for any arbitrary jet algorithm if we require that 
each factorized part be infrared hnite. If any factorized part contains infrared divergence 
for a certain jet algorithm, though the partial sum of the factorized parts is infrared safe, 
the factorization approach loses its physical meaning. 

In order to approach the problem about which jet algorithms respect the factorization 
theorem, we devise a general jet algorithm, motivated by the jet algorithm, in which we 
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introduce a parameter a. The original exclusive kx algorithm is given by 


Uij 
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(1 — cos 9 ij) 


min.{El E]) < 


( 1 ) 


where Q is the center-of-mass energy. For a given pair of partons i and j, with the energy Ei 
and Ej respectively, the parameter yij is constructed with the relative angle 6ij. Physically 
this jet algorithm requires that the relative transverse momentum of every pair of hnal-state 
partons i and j be measured. If the smallest value of yij for a given pair of partons is less 
than a resolution parameter y^, they are combined into a jet. This process is repeated until 
all pairs have yij > yc- If the relative transverse momentum is replaced by the invariant 
mass, the jet algorithm becomes the JADE algorithm. 

Here we employ the generalized exclusive jet algorithm, which is dehned as 


Vij = ^ ( 2 ) 

where the power of the energy in Eq. Q is replaced by a, and the resolution parameter y^. 
remains dimensionless. As a varies, the general jet algorithm includes all the jet algorithms 
mentioned above for specihc values of a. For example, the algorithm is obtained for 
a = 2. The JADE algorithm is similar to the case a = 1. The Cambridge/Aachen algorithm 
corresponds to a = 0, and the anti-fc-r algorithm^ is given by a = —2. We probe all the 
possible values of a and compute the jet and the soft functions to see if they are well dehned 
and free of infrared divergence. The phase space available for the collinear and the soft 
functions differs for different values of a, so does the structure of divergence. 

The next-to-leading-order contribution of each factorized part will be computed accord¬ 
ingly, in which there are at most two particles in a jet. We compute each contribution using 
the dimensional regularization with the spacetime dimension D = 4 — 2e regulating both the 
ultraviolet (UV) and the infrared (IR) divergences. Also the MS scheme is employed with 
dvr/iAg = The separation of the UV and IR divergences is important to see whether 

each factorized part is IR hnite. So far, it has been proven that the IR divergence cancels 
in many processes. And, in fact, any physical quantity should be IR safe. In this case, the 
remaining divergence is of the UV origin. Technically, we can put all the scaleless integrals 
to be zero, and also put all the poles in e as the UV poles in euv- But this is possible only 
after it is verihed that each factorized part is IR hnite. 

We intend to verify which types of the jet algorithms guarantee the IR hniteness of each 
factorized part, and which do not, using the generalized jet algorithm. As it turns out, the 
jet algorithm for a < 2 indeed yields IR-£nite jet and soft functions, hence can safely put 
all the poles as the UV poles. However, for a > 2 each jet and soft function contains IR 
divergence, and we cannot naively put all the poles as the UV poles, though the jet cross 
section is IR finite. This verihcation is possible only when we distinguish the UV and IR 
divergences, which we explicitly show here at next-to-leading order. The existence of the IR 


^ The anti-/cT algorithm in Ref. m is proposed for the inclusive jet algorithm. We refer to the exclusive 
jet algorithm with a = —2 in Eq. (l2|) as the anti-fc^ algorithm. 
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divergence in the factorized part, or the breakdown of the factorization theorem was also 
considered in Ref. [2] in a different context. 


In full QCD, only the jet cross section has been computed and the divergence is regarded 
as the UV divergence. This is justihed since the cross section is a physical quantity free 
of IR divergence. However, the factorization property in high-energy scattering of the jet 
cross section or any other physical observables such as thrust is important to disentangle 
the strong interaction and to compute high-order corrections for better accuracy. Then it 
is crucial to see if the factorization works or breaks depending on the jet algorithms. We 
address this problem by considering the IR hniteness of the jet and the soft functions in 
various jet algorithms and point out which jet algorithms allow the factorization theorem in 
the jet cross section. 


The structure of the paper is as follows: In Section the factorization of the dijet cross 
section with the jet algorithm is briefly described in the framework of SCET. In Section |ITT[ 
the phase space for the naive collinear, the zero-bin and the soft contributions is analyzed 
using the generalized jet algorithm. The structure of the phase space shows different behavior 
for different values of a, and we classify them according to the possible divergence structure. 


In Section IV the collinear and soft contributions for a > 0 are computed. These 


contributions are considered for various values of a, corresponding to the different phase 
spaces. We verify that the collinear and soft functions are infrared hnite for 0 < a < 2, while 
they contain infrared divergence for a > 2. The case a = 2, corresponding to the original 
kx algorithm is not well dehned in dimensional regularization though it can be inferred that 
each contribution contains the IR divergence from the consideration of the phase space. In 
Section |V[ the collinear and soft contributions are computed for negative a. In this case, 
there should be a jet veto needed in calculating the soft function. Each contribution is IR 
hnite for any negative a, and the generalized jet algorithm with a < 0 shows similar behavior 
compared to the cone-type and the Sterman-Weinberg jet algorithms na. In Section [VH we 


discuss the structure of divergence for various values of a. The identihcation of the UV and 
IR divergences can be inferred from the structure of the corresponding phase space. It is 
explained in detail what kind of divergence there will be a priori by examining the structure 
of the phase space. In Section VII we discuss peculiar characteristics of the kx algorithm 
with a = 2. In dimensional regularization, it is impossible to extract the poles in the jet 
and the soft functions. Furthermore they contain poles in l/(a — 2). The behavior of the 
collinear and soft parts for a = 2 is explained in detail. 


In Section [VlII[ all the contributions are added to compute the dijet cross section. The jet 
algorithm with a <2 yields hnite collinear and soft functions. For a > 2, the factorization 
breaks down, that is, each function contains IR divergence. However, if we add all the 
contributions, the IR divergence and the singularity in l/(a—2) cancel and the resultant dijet 
cross section is IR hnite. We also present the resummed results at next-to-leading logarithmic 
accuracy for the dijet cross sections with their theoretical uncertainty, and compare with 
the hxed-order results. In Section IX, we summarize the characteristics of the jet and the 
soft functions in the generalized jet algorithm and give conclusions. In Appendix, we collect 
all the ingredients for the numerical analysis. 
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II. FACTORIZATION OF THE DUET CROSS SECTION 


The factorization theorem of the dijet cross section has been presented in Ref. [13], and 
we briefly review the resnlt instead of rederiving it. The factorized dijet cross section is 
given by 

^jet , P-jjTn,©(P")(h)• (2) 

Here is the invariant-mass sqnared of the system, and cxo is the Born cross section 
for a given flavor / of the qnark-antiqnark pair with the electric charge Qj, given by 


^ 3g2 


(4) 


is the hard fnnction which is obtained from the matching of the electromagnetic 
cnrrent between the fnll QCD and SCET at leading order as 


J^ = C{Q\fi)XnYlYY^Xw, 


(5) 


and = |C'(g^,/i)p. To one loop, it is given by [T6 


H{Q^,fx) = l + 


asCp ^ oT... 

27r 




Q‘ 
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And Xn is a gange-invariant collinear qnark with a collinear Wilson line Xr, 
Y (x) is the soft Wilson line [TT] 


Ynix) = Pexp ig / dsn ■ As{sn) 


( 6 ) 

and 

(7) 


where ‘P’ denotes the path ordering. 

After redehning the collinear fields, the soft interaction is deconpled from the collinear 
interaction [S], and the factorized jet cross section in Eq. ([^ is obtained. The integrated jet 
fnnction is defined with the jet algorithm as 


Jn,©(fi) = J (ipVn,©(p^p), (8) 

where is the invariant mass sqnared of the collinear jet. The nnintegrated jet fnnction 
Jn,e is defined as 

II(O|Xnl^n)0j(A:„|Xn|O) = J ^n,© (9) 


where 0 specihes the jet algorithm. The soft fnnction 5© with the jet algorithm is given by 


Se = Y. ^'IV(0|y;F„|.Y,>e„ft(.Y,|F„tyW|0>. 


( 10 ) 


The jet fnnction J7n,© and the soft fnnction Sq are normalized to 1 at tree level, and we 
compnte them at next-to-leading order with varions valnes of a in the generalized kT algo¬ 
rithm. 
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III. JET ALGORITHMS 


We now consider the constraints imposed on the phase space by the jet algorithm in 
Eq. ([^. At next-to-leading order (NLO), there are at most two particles in a jet. The 
dominant contribution for the dijet in e+e“ annihilation is the quark-antiquark pair forming 
a back-to-back jet and a gluon is emitted. We choose the jet in the lightlike n direction, or 
in the n direction with = 0, = 0 and n -n = 2. 

In applying the jet algorithm, the virtual correction and the real gluon emission are 
included. These can be obtained by cutting the diagram in Fig. which is part of the 
matrix elements squared for the jet cross section. If we cut a single line, it corresponds to 
the virtual correction. When we cut the loop, there are two particles in the final state with 
momenta I (for a gluon) and p — I (for a quark) in the n direction. The following nontrivial 
consideration of the phase space using the jet algorithm applies to this case with a real gluon 
emission. 

The total momentum of the hnal state scales as 

= {n-p,p±,n-p) = {p-,p±,p+) ~ Q(1,A,A^), (11) 

where A is the small parameter in SCET. We choose the n direction as the direction of p 
such that p_L = 0, and p- = Q. The collinear gluon momentum scales as 

P = (L,/x,/+)~g(l,A,A2). (12) 

With this power counting, the collinear momenta of the quark and the gluon can be written 
as 

= (Q - -h,pVQ - 4), = (/-, 4,4), (13) 

with their energies 

Eq = 2 ^Q ~ + P^IQ ~ W)i Eg = -{I- + 4 )- ( 14 ) 

And the invariant-mass squared p^ is given by 


i 



FIG. 1. Particle configuration and the momentum assignment in constructing the phas space. 
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The phase space constraint from Eq. ([^ for the jet fnnction changes with the sign of a. 
In terms of the glnon momentnm I and Q, the jet algorithm for a > 0 is given as 

22-a I la-l 


y,.-> 


'> 2 — 0 : 




< Vc, 


0 < i_ < |, 

Q 


(l - g) p < 2/c, ^ < /- < Q. 


(16) 


/_ 


The nontrivial jet algorithm appears when the resolution parameter is of order in 
SCET. Hence we choose yc ~ 0{X‘^) from now on. 

By letting l_ = Qx, /+ = yy, the phase space is constrained bj|^ 


Vga 


22-a 

22-« 


Ct — 1 

/i yx 

Q{i-xy 

Q X 


< Vc, 
a-2 

- < Vc, 


0 < X < ^, 
\ < X <1. 


Solving for y, it becomes 

y < = /.(^). (0 < I < i), 

y < 2“"V^a:(l - xf-'" = (^ < X < 1), 


(17) 


(18) 


where the boundaries are defined in terms of fi{x) and f 2 {x) respectively. 

The shape of the phase space in the (/_,/+) space with a > 0 changes drastically for 
the following five cases: (i) 0 < a < 1, (ii) a = 1, (hi) 1 < a < 2, (iv) a = 2, (v) a > 2. 
The corresponding phase spaces are shown in Fig. The hgures in the hrst column are 
the phase spaces in the (/_,/+) space, and those in the second column are those in the 
(Z_,l^) space. The structure of the divergence can be inferred from the phase space in the 
(/_,/+) space. But when it becomes ambiguous, the phase space in the (/_,1^) space plays 
a complementary role in deducing the origin of the divergence [TTl [T8] . This point will be 
discussed in detail in Section |V^ after the explicit result is presented. 

For 0 < a < 1, /+ approaches zero at the endpoints Z_ = 0 and Q. At a = 1, /+ becomes 
a nonzero constant at Z_ = 0, while it approaches zero at /_ = Q. For 1 < a < 2, /+ diverges 
as Z_ approaches zero, while it becomes zero at /_ = Q. For a = 2, diverges at /_ = 0, 
while remains finite at /_ = Q. For a > 2, /+ diverges at both endpoints /_ = 0 and Q. The 
different behavior of the phase space near the endpoints gives a clue to whether the naive 
contribution to the jet function contains UV or IR divergence. 

In the hrst column of Fig. Fig. (a) and (b) corresponding to a < 1, /_ and Z+ never 
reach inhnity, and we expect that the naive collinear contribution has IR divergence only. 
In Fig. 1^ (c), (d) and (e) corresponding to a > 1, can reach inhnity and we can expect 
that there may be UV divergence as well as IR divergence when the integration over /_ and 
/+ is performed. However, this is not exactly true. 


^ The coefficient of the rescaling depends on whether the corresponding momentum reaches a finite value, 
say, Q, or infinity. For the momentum reaching infinity, the dimensionful coefficient is the renormalization 
scale y which is to be combined with the integral in applying the dimensional regularization. See the 
zero-bin contribution or the soft function below. 
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FIG. 2. The phase space in the general kx algorithm with a > 0. The first column shows the 
(Z_,Z+) phase space, while the second column shows the (/_,lj^) space, (a) 0 < a < 1, (b) a = 1, 
(c) 1 < a < 2, (d) a = 2, and (e) a > 2. The functions fi and /2 describe the boundaries of the 
phase space for 0 < < <5/2 and Ql2 < < Q respectively, as defined in Eq. (18). 


In most cases the UV and IR divergences can be expected by looking at the (Z_, Z+) phase 
space alone. But there are some occasions where it is ambiguous. Then we turn to the 
(/_,1^) phase space. For example, if we consider the hnal form of the integral in computing 
the collinear part with respect to /_ and 1^, the second hgure in Fig. (c) shows that 1^^ 
never reaches inhnity, hence there should be only IR divergence for 1 < a < 2. For a > 2, 
both hgures in Fig. (e) indicate that there are both IR and UV divergences. In the second 
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figure in Fig. (d), we expect that there is only IR divergence, but the collinear and the 
soft functions are not well-dehned for a = 2, and it needs more explanation, as we will see 
below. 

In order to avoid double counting, we subtract the soft limit from the naive collinear con¬ 
tribution, which is referred to as the zero-bin subtraction [T^]. For the zero-bin contribution, 
we take the soft limit of the gluon (/_, l±, l^) ~ Q{X^, A^, A^) and rewrite the hrst equation 



FIG. 3. The phase space for the zero-bin contributions with a > 0 in the same format as Fig. 
(a) 0 < a < 1, (b) a = 1, (c) 1 < a < 2, (d) a = 2, and (e) a > 2. 
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in Eq. (18). Then the phase space for the zero-bin contribution is given as 


y 


0 

99 


7 /““I 




< Vc, 


( 19 ) 


which can be written, by putting /_ = fix, /+ = p?/, as 



( 20 ) 


The phase space for the zero-bin contribution corresponding to each case in Fig. |^is shown 
in Fig. 1^ As can be seen clearly in Fig. the zero-bin contribution contains both UV and 
IR divergences. However, the IR divergence is cancelled in the sum of the naive collinear 
contribution and the zero-bin contribution for a < 2, while the cancellation is incomplete 
and contains IR divergence for a > 2. 

For the soft function, i/ij satishes the following relations with the appropriate power 
counting. 


y 


S 

Q9 


,2-a4(4 + ^-)" ^ 
,2- 4-(4 + /-)“-^ 


< yc, 

< yc, 


n jet, 
n jet. 


( 21 ) 


By replacing /_ = /ix, /+ = /i?/, the boundary from the n jet is given by 


y = 9{x) = -x+(^Y ^ 

oc 

where a is a dimensionless parameter given as 



( 22 ) 


(23) 


The boundary from the n jet is given hy y = g~^{x). The phase space for the soft function is 
shown in Fig. As can be seen in Fig. since the region near the origin is totally covered, 
no jet veto is necessary for a > 0. However, the jet veto is needed for a < 0, as in the 
Sterman-Weinberg jet algorithm |15j . 

For a < 0, the minimum in Eq. (|^ is switched compared to the case a > 0. It means 
that the boundary functions are also reversed. That is, the constraint of the phase space for 
the jet function is given by 

y < 2 “" 4 c- 2:(1 - xf-"^ = / 2 (x), (0 < x < ^), 

jJj Zi 

y < 2“"4c^2:^""(1 - xf = /i(x), < x < 1). (24) 

The phase space for the zero-bin contribution is given by 


7/ < 2" ^y^x, 


(25) 
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0 < a < 1 



FIG. 4. The phase space for the soft function with a > 0 corresponding to Fig. The jet veto 
is not necessary for a > 0, and Sa = (a) 0 < a < 1, (b) a = 1, (c) 1 < a < 2, (d) 

a = 2, and (e) a > 2. 


and that for the soft function is given by 

y < n jet, 

1 

V > -— X, n iet, 

X + y < 2/3—, jet veto. (26) 

/i 
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FIG. 5. The phase space in general kT algorithm with a < 0. (a) naive collinear contribution (b) 
the zero-bin contribution (c) the soft function. 


In contrast to the case with a > 0, the jet veto is needed here since the phase space near 
the origin is not completely covered by the jet algorithm alone. And we take (3 ~ 0{X^) for 
dehnite power counting. 

The phase space for the naive collinear part, the zero-bin contribution, and the soft 
function for a < 0 is shown in Fig. The shapes of the phase spaces are similar to the 
Sterman-Weinberg algorithm or the cone-type algorithm [TB]. And we expect that the jet 
and the soft functions are IR hnite for all negative values of a. 


IV. GENERALIZED kr JET ALGORITHM WITH a > 0 

A. Jet function 

The Feynman diagrams for the jet function in the n direction, are shown in Fig. The 
dashed line is the cut, therefore Fig. (a) is the virtual correction, while Fig. (b) and 
(c) are the real contributions. The mirror images of (a) and (b) are omitted in Fig. but 
they are included in the computation. Each diagram is accompanied by the corresponding 
zero-bin contribution, which should be subtracted to obtain the collinear jet function. 



FIG. 6. Feynman diagrams for the jet function at one loop (a) virtual correction (b) real gluon 
emission from the Wilson line (c) real gluon emission. 
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The dimensional regularization states that 


POO 1 1 

-, (27) 

“'0 euv ^iR 

where I is a momentum variable. The integral vanishes identically if we do not distinguish 
the UV and IR poles since the integral is a scaleless integral. However, we distinguish the 
UV and IR poles since we are interested in identifying the sources of the divergence. 

The virtual correction in Fig. (a) is independent of the jet algorithms and the net 
collinear contribution Ma is given by the naive collinear contribution Mq, subtracted by the 
zero-bin contribution M°: 



1. The case 0 < a < 1 


The naive collinear contribution for 0 < a < 1 in Fig. (b) is given 

agCp //i 


as 


M,= 


ri/2 _ ./i(x) _ 

2. r(i - <) yqi vh i ■■ 

/*! /*/2(^) 

/ — x) / dyy~^~‘^ 

Ji /2 Jo 


OisCp 


2ti Lq, 




2 " efu eiR 
2 


In 




)a—2 


VcQ^ 

1 /i^ /A 


1 , 2 

2 2“-2|/,g2 


a In 2 


+ 


TT 

12 


(29) 


This explicit calculation supports the expectation from the phase space in Fig. (a) that 
the poles are of the IR origin. This can be also seen by the fact that the integral of Mb in 


Eq. (29) exists only when e < 0. 


The corresponding zero-bin contribution is given by 






s'-'F 


a 7 Be 






dyy 


-i-t 


(30) 


where fi{x) is the zero-bin limit of fi{x) and is given in Eq. (20). From Fig. since 
/+ can both approach zero and infinity, the zero-bin contribution contains both UV and IR 
divergences. In order to separate the UV and IR contributions, the integral is decomposed 
to be written as 


= 


a,C 




3 7ec 


/?U) 


2 . r(i-.)U„ 
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+ / dxx ^ 

Jri 

ctgCp / 1 


dyy 


-l-e 




dyy 


-l-e 


( -)[-( 


H-h In ■ 




+ — + In - 

27r '^euv cir^ Lq, _ 2 2"“2j^^Q2y qJ’ 


(31) 


where 77 is a small positive number to separate the IR and UV regions. The hnal result is 
independent of y. This decomposition looks complicated, but it is arranged in such a way 
that the pole from the integral including the inhnity (zero) is of the UV (IR) origin with 
e > 0 (e < 0) |I5]. 

And the net collinear contribution from Fig. (b) is given by 


M, = M,- 


(32) 


dgCp [ 1 / 1 1 

- 9 V 




In 

1 . /I 1 N / 1 




-^In' 


27r La-2V euv2“-VQ' 2'“ 2“-^^" ' 12^ 


+ 19) 


+ ln^) + (4-a)(l-^) -aln2 


eiR 2" "^ycQ^ '"euv em^'^euv Q 

The naive collinear contribution from Fig. (c) is given by 


71 

12' 




277 F(1 — e) '^Q 




dyy 


-l-e 


ttgCp 


-In 




277 V 2eiR 2 2 “- 27 /,Q 2 


+ 


a — 3 OL 


+ -ln2 , 
2 2 ' 


(33) 


while the zero-bin contribution is subleading and suppressed. Adding all the terms with the 
wavefunction renormalization and the residue at one loop 


+ “.Cf 


(— - !-). 

V9e:TT.. 


277 ^2euv 2eiRU 

the collinear contribution for 0 < a < 1 at order as is IR finite and is given as 


(34) 


A-C,”*' = 2(M„ + M,) + M, + Zf > + Rf 


OisCp 




In 


/i 


-hU 


y 


TT^ 
+ ^ 


27r La-2V e^y epy 2»-^y,Q^ 2--^y,Q^ 6 

,2 . 1 Q _ 


In 




2 Veuv 2“-2i/,g 


+ 


-In 2 -I- (a — 4) 

2 ^ ^ 6 J 


(35) 


By removing the UV poles, the jet function for 0 < a < 1 at order is given by 


J!{^\Q^yc,y) = 


asCp 


277 Lq, — 2 


— f-ln^ 

- 9 V 


y 


77 ^ \ 3 

+ ^ ) + In 
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2“-27/^g2 g J 2 2“-27/^g2 


77^ 


13 — 3 q; 3a , , 

-n2+ a-4 — 

2 2 ^ ^ 6 J 


(36) 
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2. The case a = 1 


Note that the phase space for a = 1 (naive and zero-bin) allows the same method of 
integration with the same specification of the poles as in the case 0 < a < 1. The only 
difference is that /+ approaches a nonzero, finite valne as /_ approaches zero in the phase 
space for the naive collinear contribntion, hence the case a = 1 can be inclnded in Eq. (35). 
The case a = 1 is very similar to the JADE jet algorithm as far as the phase space is 
concerned. 

In the JADE algorithm, the invariant mass of every pair of the final-state partons i 
and k should be less than where j is some fraction (of order in the context of SCET) 
to form a jet. At next-to-leading order, the JADE algorithm is written as 

/+ < j{Q — 1-), collinear, 

/+ < jQ, zero-bin, 

/_ < jQ, /+ < jQ, soft. (37) 

If we compare this phase space with that of Fig. (b).| (b) and (b) for a = 1, the 
basic structure near the endpoints of the phase space is the same, and the only difference 
lies in the middle region of Fig. (b). Therefore we expect that the divergence structure for 
the case a = 1 and the JADE algorithm is the same except the finite part in the dijet cross 
section. In fact, performing the calculation using the JADE algorithm, the collinear part is 
given as 
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which is consistent with the result in Ref. m- And the corresponding collinear part for 
a = 1 from Eq. (35) is given by 
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(39) 


If we put j = i/c/2, all the UV divergence and the logarithmic terms are reproduced except 
the finite terms, which is why we say that the generalized jet algorithm with a = 1 is 
basically identical to the JADE algorithm. 

From Eqs. (38) and (39), the jet functions in the JADE algorithm and in the case with 
a = 1 are given as 
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(40) 


3. The case 1 < a < 2 

By looking at the phase space for the naive collinear contribution in Fig. (c) in the 
(/_,/+) space, one might naively expect that there could be a mixed UV and IR divergence 
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since /+ reaches infinity as /_ approaches zero. However it becomes clear when we look at 
the phase space in the (/_,1^) space. It means that, when we perform the integral over I, 
the is first performed and the remaining integral is to be performed with respect to /_ 
and [HI [18]. Then the divergent behavior is determined by 1^^ = If the divergence 
occurs when I^ —>■ cx), it is the UV divergence. On the other hand, if the divergence occurs 
when I^ —0, it is the IR divergence. The constraint for I^ from the jet algorithm is given 
by 

li = 4 L<<^ V 2 ’ 

y 2°‘~‘^Qlj,ycx'^{l — \ < x <1. 

which approaches zero (1 < a < 2) as x —)■ 0. It can be also seen in Fig. (c) with the plot 
in the (/_,1^) space, 1^ never approaches infinity. Therefore the poles should be of the IR 
origin. 

The collinear contributions computed for the case 0 < a < 1 can be exactly performed 
with the same prescription on the signs of e, hence giving the same result for 1 < a < 2. It 
also turns out that 

hm = Inn Ml<r^. (42) 


Therefore Eq. (35) can be extended to the case 1 < a < 2, and it is continuous at a = 1. 


4- The case a > 2 


The situation is completely different for a > 2. In the naive collinear and the zero-bin 
contributions, the phase space shows that 1^ can reach zero or infinity. Therefore we expect 
that the collinear contributions contain both UV and IR poles, which makes the factorization 
obsolete because the jet function and the soft function are not IR finite. However, as will 
be seen later, the sum of the jet functions and the soft function is IR finite. 

In the real gluon emission, the expression for is given by 
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But the function /i(x) describing the boundary of the phase space for 0 < /_ < Q/2 diverges 
as Z_ approaches zero by noting that the function /i(x) is proportional to x^“"(l — x)^. It 
requires the modification of the first integral in Eq. (43) as 
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This manipulation is necessary since the remaining integral converges for e > 0, hence the 
integral over y should produce euv- The boundary of the phase space, f 2 {x) for Q/2 <l_ < 
Q is proportional to a;(l — and it also diverges as Z_ approaches Q (x —)■ 1). However, 

the integrand is proportional to 1 — x, and there is no UV divergence involved. The naive 
collinear contribution is given as 


~ _ asCp r 1 / 1_1, 2 

27r ely euv 2 ^ ^ 12) 

1 1 1/1 

H —2 -1-^ (-) In — + In — 

Hr cuvHr Hr ^cir euv^ Q 2 ^ycQ^ 

2 2 

+ 2 - 21n2 - ^ - (a - 2)(l + ln2 - . 

D ^ 12^J 


(45) 


In the zero-bin contribution, a similar modihcation of the integral is necessary. It is given 


as 
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The net collinear contribution is given by 
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Fig.g(c ) yields the same results as in the case 0 < a < 2. Summing all the contributions, 
the net collinear part is given as 
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In this case, the collinear part is both UV and IR divergent. Therefore it does not have 
any physical meaning for a > 2 . Compared with the collinear part for 0 < a < 2 given in 
Eq. (35) , if we put euv = ^ir = the result is the same. However, we stress that this is 
misleading because the collinear part for 0 < a < 2 is IR finite and there are only UV poles, 
while that for a > 2 contains both UV and IR divergences. In this case, simply putting 
euv = eiR = e ruins the physical implication on the divergences. But, as we will see later, 
the sum of the IR divergent collinear part and the soft part, which is also IR divergent, is 
IR finite. Therefore the combined collinear and soft part is physically well-defined, but the 
factorization is not realized in a strict sense. 


5. The case a = 2 

The jet algorithm with a = 2 corresponds to the exclusive /c-r algorithm. A serious 
trouble in this case is that the jet function has a pole at a = 2 and is not well-defined. 
This is to be contrasted with the case a = 1. For a = 1, the structure of the phase space 
remains continuous for 0 < a < 1 and 1 < a < 2 , hence the same divergence structure 
follows. However, for a = 2, the shape of the phase space abruptly changes from 0 < a < 1 
to a > 2. The structure of the divergence also abruptly changes accordingly. There is only 
IR divergence for 0 < a < 1, while there are UV, IR divergences for a > 2. Therefore the 
jet function for a = 2 cannot be obtained by taking the limiting case as a —)■ 2 from either 
side. This is manifested by the existence of the pole in l/(a — 2). 

However, as we will see later, the dijet cross section in the kj- algorithm is finite, but 
we emphasize that the factorization breaks down for a = 2. Since the behavior of the jet 
function, as well as the soft function, is peculiar, it is worth mentioning what happens at 
a = 2 in detail later. 

B. Soft function 

The Feynman diagrams for the soft function at one loop is shown in Fig. in which the 
hermitian conjugate is omitted. The phase space for the soft function is shown in Fig. 



(a) (b) 


FIG. 7. Feynman diagrams for the soft function at one loop (a) virtual corrections (b) real gluon 
emission. 
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The prescription for the UV and IR divergences is the same as in the case of the jet function 
according to the possible values 1^. 

The virtual correction in Fig. (a) is given as 
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(49) 


which is independent of the jet algorithm. We compute the real gluon emission from Fig. 
(b) for various values of a > 0. 


1. The case 0 < a < 1 


The structure of the phase space for the soft function is shown in Fig. (a), and it is 
apparent that both the UV and IR divergences may exist. The contribution of the real gluon 
emission is given by 
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where A = (2^ is the point where two boundaries meet, and g(x) describes the 


boundary for the n jet, and is given in Eq. (22). We divide the phase space symmetrically 


with respect to the line y = x, and compute only the upper half, hence the factor 2 is 
multiplied. By carefully separating the UV and IR divergences, the amplitude for the soft 
function, including the hermitian conjuagates, is given by 

= 2(^, + S,) (51) 
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where a = 2°‘~‘^yc{Q/ ij)°‘■ The soft function at one loop for 0 < a < 1 is given as 
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27r La(a — 2) 

The IR divergence from the virtual correction is cancelled by that in the real gluon 
emission, and there remains only UV divergence. This fact can be seen transparently from 
the argument of the phase space |T8]. The virtual contribution corresponds to the same 
integral as for the real gluon emission, but the integration region covers all the phase space, 
with a minus sign. Therefore the total soft contribution is given by the integral over the 
unshaded phase space in Fig. |^with a minus sign. Since the unshaded region never touches 
the IR region, the soft contribution contains only UV divergence. 


2. The case a = 1 


The calculation procedure is the same for a = 1, hence the result can be extended to the 
case a = 1. The soft part at a = 1 is given as 
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while the soft part in the JADE algorithm is given by 
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( 54 ) 


which is identical to the soft part with a = 1 by replacing j = yd^^,. This should be true 
because the jet algorithm for a = 1 in the soft part is identical to the JADE algorithm with 
j = ycd- The soft functions in these two algorithms at next-to-leading order are given as 
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3. The case 1 < a < 2 


The contribution from the real gluon emission is given as 
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(56) 


It turns out that the result is the same as that in the case with 0 < a < 1. 


4- The case a > 2 


The left-hand-side plot in Fig. (e), Z+ approaches inhnity as Z_ goes to zero. Therefore 
we may expect that there are both UV and IR divergences. However, the plot of 1^^ with 
respect to /_ in the right-hand side shows that 1^ can never reach infinity, hence the real 
gluon emission involves only the IR divergence. In order to handle the divergence correctly, 
we write Sb as 
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The amplitude for the soft function is given by 
AC?J = 2(S„ + SO 
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Note the difference of the soft amplitudes between the cases with 1 < a <2 and a > 2. 
If we naively put euv = ^ir = e, the result is the same. However, if we distinguish the UV 
and IR divergences, the meaning of the results is totally different. The case with 1 < a < 2 
contains only the UV poles, while the case with a > 2 is not IR hnite. Therefore, when 
a > 2, the soft function is not physically meaningful. However, if we add the soft amplitude 
with the corresponding collinear amplitude, Eq. (48), the IR divergence cancels out. It 
means that the dijet cross section can be computed with IR hniteness, but the factorization 
breaks down for a > 2. 

One more comment is that the results are also divergent at a = 2. The soft function 
for a = 2 is not well dehned in dimensional regularization. However, if we add the soft 
amplitude with the corresponding collinear amplitude, the divergent term in a, proportional 


to l/(a — 2) also cancels out. This will be discussed in detail in Section VH 


V. GENERALIZED kr JET ALGORITHM WITH a < 0 

The generalized algorithm for a < 0 is similar to the generic cone type jet algorithm 
or the Sterman-Weinberg jet algorithm as far as the shape of the phase space is concerned 
[15]. The phase space for the naive collinear contribution, the zero-bin contribution and the 
soft function is shown in Fig. We can extend the jet algorithm to a = 0 since the phase 
space for a = 0 is obtained by taking the limit a —)■ 0 from the case with a < 0 because 
a jet veto is needed. The jet algorithm with a = 0 is known as the Cambridge/Aachen jet 
algorithm. 


A. Jet function 


The naive collinear contribution from Fig. (b) is given 
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(59) 


For the zero-bin contribution, the boundary for the phase space is obtained by taking the 
zero-bin limit of f 2 {x), which is given by 

f^{x) = 2^-^y^x. (60) 


The zero-bin contribution is given as 
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From the real gluon emission in Fig. (c), the naive collinear contribution is given by 
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Compared to the corresponding calculation for a > 0 in Eq. (33), the only difference is in 
the coefficient of In 2. 

Combining the wave function renormalization and the residue, the collinear amplitude 
for a < 0 is given as 
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And the jet function at one loop for a < 0 is given as 
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B. Soft function 


The contribution from the real gluon emission is given by 
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With the virtual correction, the amplitude for the soft function is written as 
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and the soft function at one loop for a < 0 is given by 
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In the limit a —>■ 0, which corresponds to the Cambridge/Aachen jet algorithm, the soft 
fnnction is given as 
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( 68 ) 


Note that the limit a —)■ 0 should be taken from a < 0. For a > 0, there is a pole in 1/a, 
and it may be dangerous to take the limit a —)■ 0. However, when the limit a —)■ 0 is taken, 
it coincides with the phase space for a < 0 and the jet veto is needed. Therefore the soft 
function at a = 0 should be taken from the soft part with a < 0. 


VI. PHASE SPACE AND STRUCTURE OF DIVERGENCE 


The structure of divergence can be inferred from the phase space argument [HI |T8] . Since 
the phase spaces for different a yield different shapes of the phase space, we will discuss and 
compare the cases for various values of a. 

Let us hrst consider the phase space for the naive collinear contribution in Fig. In 
the hrst column, the phase space is plotted in (/_,/+) space, and the second column in 
(/_,1^) space. When we compute the naive collinear contribution, there is an integration 
with respect to the gluon momentum 1. If we hrst integrate over 1^, which takes care of the 
delta function in the integrand, the remaining integral is to be performed with respect to /_ 
and /+, where the phase space in the hrst column is relevant. If we hrst integrate over 
the remaining integral is to be performed with respect to Z_ and 1 _l, where the phase space 
in the second column is appropriate. 

In (/_,/+) space, the UV (IR) divergence comes from the region where /_ and /+ ap¬ 
proaches inhnity (zero). If one of them approaches inhnity, while the other approaches zero, 
it is possible that there is a mixed UV-IR divergence. We can infer the structure of diver¬ 
gence in (Z_, Z+) space from the hrst column of Fig. For a < 1, it is clear that, if there is 
divergence, it should be of the IR origin from Fig. (a) and (b) since Z_ and Z_|_ can never 
reach inhnity. It is explicitly verihed in Eqs. (29), ( [M| ). This is also evident from the second 
column of the corresponding phase spaces since Z_ and l^ never reach inhnity. With the 
same argument, when a > 2 in Fig. (e), both plots imply that there are UV, IR and mixed 
divergences. 

However, the two plots lead to seemingly inconsistent results in Fig. (c), which corre¬ 
sponds to 1 < a < 2. In (/_,/+) space, /+ can reach inhnity while /_ cannot. Therefore 
we naively expect that there might be IR and mixed divergences. On the other hand, in 
(/_, IH space, l^^ never reaches inhnity and there should be only IR divergence. In this case, 
we have to rely on the phase space in (Z_,lH space because there is an additional relation 
Z_i_Z_ = l^^ for the real gluon emission. 


From the phase space constraint in Eq. (16), it reads for /_ < Q/2 
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and indeed /+ can approach infinity as /_ approaches zero for 1 < a < 2. 
bonnded by 
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which cannot reach infinity as /_ approaches zero. Since the divergence behavior is deter¬ 
mined by the behavior in in 77 — 2 dimensions, and = /+/_, the inference on the 
divergence behavior is more reliable in (/_,1^) space rather than in (/_,/+) space. This 
can also be confirmed in Ref. IZH, in which the jet fnnction has been compnted with the 
transverse momentnm fixed. Since the transverse momentnm is fixed at a finite valne, there 
is no UV divergence, bnt IR divergence. The compntation was performed for the inclnsive 
jet fnnction with the transverse momentnm fixed, hence there is no constraint on the phase 
space as in the jet algorithm. However, the resnlt snpports that there is only IR divergence 
when 1^ never reaches infinity. As can be seen in Fig. (a) to (c), the phase spaces for 
1 < a < 2 in (/_, 1^) space shows similar shapes and there is only IR divergence. 

For the phase space of the zero-bin contribntion in Fig. all the plots show that the 
zero-bin contribntion should contain both UV and IR divergences. For the soft contribution 
shown in Fig. we can conclude that there are both UV and IR divergences for 0 < a < 2, 
while there is only IR divergence for a > 2, as can be seen in Eq. (57). However, we can 
apply the phase space method as in Ref. 


The virtual correction for the soft part has the 
same integral as in the real gluon emission except the minus sign, and it covers the entire 
phase space. Therefore the soft part is obtained by integrating the matrix element over the 
unshaded phase space in Fig. |^with a minus sign. Then the soft part for 0 < a < 2 has only 
UV divergence since the unshaded region never reaches zero. It means that the cancellation 
of the IR divergence occurs between the virtual and real contributions. On the other hand, 
for a > 2, we expect that there is not only UV divergence, but also IR divergence, which 
breaks the factorization. 

The case with a = 2 (the kr algorithm) is tricky because it contains the pole in l/(a —2), 
and it cannot be computed in dimensional regularization. However, from the phase space 
diagrams, we expect that the naive collinear contribution contains only IR divergence, while 
the soft contribution contains both UV and IR divergences. This will be discussed in detail 
in the next section. 


VII. COMMENT ON THE kr ALGORITHM WITH a = 2 


We have considered the jet and the soft functions in the generalized jet algorithm by 
varying a. It turns out that the jet algorithm for a < 2 respects the factorization, that is, 
the jet and the soft functions are IR finite, and it breaks the factorization for a > 2 since 
each part contains IR divergence. The case with a = 2 shows a peculiar behavior and it is 
worth mentioning a few characteristics. 

At a = 2, in the kr algorithm, there is an uncontrollable integral of the form 



( 71 ) 
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in the naive collinear, the zero-bin and the soft contributions, even though the spacetime 
dimension is extended to D = 4 — 2e. It is also noted in Ref. (2Uj . 

If we use the offshellness to regulate the IR divergence in the soft Wilson line by replacing 
the propagator, the soft part at a = 2 is given by 
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where Aj = pf/Q represent the offshellness of the quark and the antiquark. And this is 
consistent with the result in Ref. [20]. The result clearly shows that the soft part contains 
the UV and IR divergences. If we look at the plots for the phase space in Figs. |(d). 
1^ (d) and 1^ (d), it can be inferred that the naive collinear contribution contains only IR 
divergence, the zero-bin and the soft contribution have both UV and IR divergences. 

In spite of the possible argument on the existence of the UV and IR divergence in the jet 
and the soft functions for a = 2, we interpret the result differently in the sense that they 
are not well dehned for a = 2. If we consider the jet and soft functions as a function of a in 
the generalized kx algorithm, they are not continuous at a = 2 in the sense that 


lim (M”f, lim (M“f, (73) 

0 .^ 2 - a^2+ 

When we look at the plots for the phase space, for example. Fig. the phase space allows 
hnite 1^ only for a < 2, but 1^ can reach inhnity for a > 2. Therefore the divergence 
structure abruptly changes at a = 2. This discontinuity appears as a pole in a — 2 in our 
treatment of the generalized jet algorithm as a function of a. In contrast, the behavior of the 
jet function and the soft function near a = 1 is smooth. The shape of the phase space in the 
(/_,/+) space changes near a = 1, that is, /+ approaches a hnite value instead of zero near 
/_ = 0. However, in the (/_,1^) space, there is no abrupt change in shape. Therefore the 
jet algorithm is continuous at a = 1, hence the factorization works in the JADE algorithm, 
too. 

Interestingly enough, the pole in a — 2 as well as the IR divergence cancels in the sum 
of the collinear part and the soft part to yield the jet cross section. Therefore the jet cross 
section is IR hnite even at a = 2. We may conclude that the kx algorithm works hue, but it 
is misleading. The main point is that the factorization, which lies behind the jet algorithm, 
breaks down for a>2 including the kx algorithm. 

However, note that the claim that the kx algorithm breaks factorization may be too 
strong because the separation of the UV and IR divergence is spoiled by the existence of the 
pole in a —2. In dimensional regularization, the extraction of the divergence is impossible at 
a = 2, but there may be other ways to obtain the divergence unambiguously. For example, 
in Ref. [22|, the authors have tried to regulate this integral with the rapidity renormalization 
technique. It may turn out to be useful, but we will not pursue in this direction here, and 
our claim stands as it is as far as we use the dimensional regularization to extract UV and 
IR divergences. 
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VIII. RENORMALIZATION GROUP BEHAVIOR OF THE DUET CROSS SEC¬ 
TIONS 


Let us collect all the results for the resummation of the large logarithms of the form 
a” In^ Dc with n < k < 2n through the renormalization group behavior. At one loop, the 
hard function is given by 






(74) 


The jet function and the soft function for 0 < a < 2 at one loop are given as 
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The jet function in the n direction is the same as ■ The dijet cross section in Eq. (3) 
at next-to-leading order (NLO) for 0 < a < 2 is written as 


(1) _ OtsCp 

^0<a<2 — ^0 2^ 


■ - hT ( 2 "- 2 ^j _ 3 ln(2^-%) + 5 - 3a - 3a In 2 + 
a 


3a 


-TT 


(76) 


For example, the dijet cross section for a = 1, which is close to the JADE algorithm, is 
given by 
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while the cross section exactly in the JADE algorithm is given by 
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This agrees with the full QCD calculation 

For a < 0, the jet function and the soft function at one loop are given as 
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The dijet cross section at NLO for a < 0 is given as 
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An interesting case is the result of the anti-/cr jet algorithm with a = —2. The dijet cross 
section in the anti-Zc^ algorithm at order o;* is explicitly given by 
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Also the Cambridge/Aachen jet algorithm corresponds to the case a = 0, and the dijet cross 
section in this jet algorithm is given as 
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For a > 2, the factorization breaks down since the jet and the soft functions contain IR 
divergence. However, when the jet function and the soft function are added, these divergent 
terms cancel and the hnite dijet cross section is obtained. The sum of the collinear and the 
soft parts is given as 
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Note that only this combination is free of IR divergence as well as of the pole in a — 2. 
Therefore the dijet cross section has a hnite limit as a —)■ 2. The dijet cross section for 
q; > 2 is given as 
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Note that the dijet cross section is the same for a < 2 and a > 2, and it is continuous at 
a = 2. Therefore, for the kx algorithm, the dijet cross section at order as is given as 
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The hxed-order dijet cross sections for various values of a are independent of the renor¬ 
malization scale /i. However, they contain logarithms of yc- When y^ ~ 0{X‘^) becomes 
small, the hxed-order results are not to be trusted and the resummation of the large log¬ 
arithms should be performed. This is achieved by solving the renormalization group (RG) 
equations for the hard, jet and soft functions. We present the solution of the RG equation 
and show the numerical results of the resummation ehect with the theoretical uncertainties 
at next-to-leading logarithm (NLL) order. 

The RG equation for the hard, jet and soft functions is of the form 
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where y) is the corresponding function and Ui = Q, yj, ys respectively for i = H, J, S. 
The solution of the RG equation can be expressed in terms of the following quantities which 


27 



are defined as 
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Then fi) is written as 
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where the scales /Xj for i = H, J, S are the hard, jet and soft scales respectively 

The (3 function and the anomalous dimensions can be expanded as a series in ag as 
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The anomalous dimension for the hard coefficient at order is given by 
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The anomalous dimensions of the jet and the soft functions for 0 < a < 2 at order ag are 
given as 
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where the jet scale fij and the soft scale fis are dehned as 
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For a < 0, they are given as 
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It can be explicitly verihed that '^h + 27 j + 75 = 0 for a < 2, and the dijet cross section 
is independent of the renormalization scale. However, the anomalous dimensions for the jet 
and the soft functions for a > 2 are not dehned since the corresponding amplitudes include 
IR divergence. In this case the anomalous dimension of the combined collinear and soft 
parts is given by 

agCp^ ,, 


729+s = ^(-41n^+6). 
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which exactly cancels that of the hard function 

The anomalous dimensions for 0 < a < 2 can be cast into the form 
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where rcusp(Q;s) is the cusp anomalous dimension of the hard function H{Q, fi). This form is 
obtained since Jh + ‘^'yj + 13 = 0 , and since the cusp anomalous dimensions are independent 
of a, that is, the jet algorithm. To NLL order with large logarithms, we need the cusp 
anomalous dimension to two-loop order, the remaining anomalous dimensions to one loop 
order and the hard, jet and soft functions at tree level m The details for the numerical 
analysis are provided in Appendix. 

The hard, jet and soft scales can be inferred by requiring that the logarithmic terms 
in the fixed-order results are not significant. The hard scale is set by = Q, and it is 
varied from /iH/2 to 2fiH to evaluate the theoretical uncertainty. The jet scale is given by 
fxj = for a < 2, and is also varied from /ij/2 to 2/ij. 

The choice of the soft scale is more delicate. For 0 < a < 2, from the argument of 
the logarithm in the soft function in Eq. (75), the soft scale is set as fis- From the power 
counting, we require that /ij ~ XQ and ~ where A is a small parameter in SCET. 
For a < 0 , the soft scale is set by = ‘2[3Q. As discussed in the jet algorithm, (3 is also 


of order A^. Since A^ is set by the jet scale as 2 " , we require that 2/3 


')a—2 


i/c in the 


numerical analysis. We also vary the soft scale from p. 5/2 to 2fis to evaluate the theoretical 
uncertainty. 

Fig. 1^ shows the theoretical uncertainty of the jet cross sections for three values of a at 
NLL as we vary the jet and soft scales. The cross sections are normalized by dividing the 
jet cross section by the total cross section at order which is given by 


aS = ao(l + ^^). (97) 

Therefore the ratio o'jet/o'tol represents the fraction for the two jets. The values of Q and 
He are chosen as 1000 GeV, and 0.1 respectively for illustrative purposes. The hard scale is 
hxed at Q. And the cases with a = 1, 0 and —2 are shown, which correspond to the JADE, 
the Cambrige/Aachen, and the anti-fcr algorithms respectively. 

The upper plots show the dependence on the jet scale /ijet as it varies from /ij/2 to 2/ij 
while the soft scale /Xsoft is hxed at /xg. The lower plots show the dependence on the soft 
scale /isoft as it varies from /i 5/2 to 2 / 1 ^ while the jet scale /Xjet is hxed at /ij. As can be seen 
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FIG. 8. The dependence of the resummed dijet cross sections on the renormalization scales for (a) 
a = 0 (b) a = —1 (c) a = —2. We set Q = 1000 GeV, and Uc = 0.1, and for the soft function 
p = 0 . 1 . 


in the plots, the dependence on the jet scale is mild, but the cross section is rather sensitive 
to the choice of the soft scale for a = 0 and —2. 

Fig. shows the theoretical uncertainty of the resummed and the hxed-order dijet cross 
sections. In resummed results, the hard, jet and soft scales are varied as Q/2 < /ihard < 2(5, 
/ij/2 < /ijet < 2/ij and Hs/‘^ < hsoft < 2/i5 respectively. In the hxed-order result, we set 
P-hard = Pjet = Psoft and phard is Varied from Q/2 to 2Q. The small parameter A corresponds 



FIG. 9. The dijet cross sections with (a) a = 1 (b) a = 0 (c) a = —2 corresponding to Fig. 
The bands show the theoretical uncertainties. The solid line represents the resummed cross section 
with /ihard = Q, Pjet = PJj Psoft = PS- The dashed line is the fixed-order cross section with ^ = Q. 
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to 2“ ^i/c, and the yc values on the horizontal axis indicate that we set A < 0.3. 

The hxed-order results diverge as yc approaches zero, but the resummed results are sup¬ 
pressed for small yc and converge to zero. The theoretical uncertainty at NLL order is not 
signihcantly reduced compared to the NLO hxed-order results. We expect that the uncer¬ 
tainty will be reduced at higher orders. In addition to the suppression of the resummed 
result for small y^ the hxed-order results tend to be more negative as a decreases, while 
the jet cross sections have meaning only when they are positive. As Figure shows, the 
NLO result and the NLL result are consistent except for small values of yc for a = 1. As 
a decreases, the hxed-order result is not reliable. Only after the resummation is performed, 
the jet cross sections become positive, including the behavior for small yc- It turns out that 
the cone-type algorithm, which is equivalent to the inclusive kx algorithm at NLO, and the 
Sterman-Weinberg algorithm show similar shapes for the phase spaces [TB], therefore we 
expect that the dijet cross section should be resummed in the cone-type algorithm to yield 
a meaningful answer. 

Since the dijet cross section is IR hnite in spite of the fact that both the jet and soft 
functions are IR divergent, one might expect that the resummed result for the kx algorithm 
{a = 2) can be obtained by taking the limit a —)■ 2 from the result with 0 < a < 2. 
Indeed we can obtain the resummed result and verify that the cross section approaches zero 
for small yc- However, the theoretical uncertainties blow up in this limit due to the large 
uncertainties in the jet and soft functions as a —?• 2, and there is no predictive power in the 
resummed result for a = 2- 


IX. CONCLUSION 

The generalized jet algorithm enables us to see the rich structure of the jet algorithm from 
the viewpoint of the factorization. It includes familiar jet algorithms with specihc values of 
a- The advantage of the generalized jet algorithm is to see the divergence structure of the 
jet and soft functions in a unihed way. We can see systematically how the structure of the 
divergence changes as a varies. 

The structure of the divergence can be inferred from the observation of the phase spaces 
in (/_,/_!_) space. When it is ambiguous for 1 < a < 2, the phase space in (/_,1^) space 
is helpful in identifying the divergence. In explicit computation, we have verihed that the 
expectation from the phase space is correct. 

To summarize, the jet algorithm for a <2 renders the jet and soft functions IR hnite, thus 
guarantees the factorization of the dijet cross section. Among the known jet algorithms, this 
includes the JADE algorithm (corresponding to a = 1), the Cambridge/Aachen algorithm 
(a = 0) and the Sterman-Weinberg and the cone-type algorithm {a < 0) and the anti- 
kx algorithm {a = —2). On the other hand, the factorization breaks down with the jet 
algorithms with a > 2 because the jet and soft functions contain infrared divergence. Though 
the IR divergence is cancelled in the dijet cross section and gives a hnite result, it breaks 
the factorization and the kx algorithm belongs to this category. 

The main point of the paper is that the structure of divergence in the collinear and soft 
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parts depends on how the jets are dehned. The important criterion is that each factorized 
part should be IR hnite in order to respect the factorization. And according to our analysis 
of the generalized jet algorithm, only the jet algorithm with a < 2 factorizes the dijet cross 
section in a strict sense. 

Note that we have considered the so-called exclusive generalized algorithms ra, which 
is based on the e’''e“ annihilation. The SCET version is presented in Ref. EOI- The inclusive 
recombination algorithm m is more suitable in hadronic collisions. But it turns out that 
the inclusive jet algorithm reduces to the cone-type algorithm [23] at next-to-leading order 
in which there are three hnal-state particles. In contrast to e’''e“ annihilation, the partonic 
center of energy is boosted along the beam direction compared to the hadron-hadron center 
of energy in hadronic collisions. Therefore the physical observables in hadron-hadron scat¬ 
tering, as in LHC, should be described in terms of the boost-invariant quantities such as 
the transverse momentum, the rapidity and the azimuthal angle. It will be interesting to 
probe the jet properties arising from the jet algorithms in other types of scattering. And the 
analysis of the divergence structure in various jet algorithms, which we have probed here, 
can be extended to deep inelastic scattering and hadron-hadron collisions with the boost 
invariance along the beam direction to see if the structure of the divergence is sustained or 
if the kinematics of the scattering changes it using the exclusive or inclusive generalized fcr 
algorithm. 
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INGREDIENTS IN OBTAINING THE NLL RESULTS 


The QGD /3 function and the cusp anomlous dimension in the MS scheme are given as 
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The coefficients of the anomalous dimensions in the hard function are given by 
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For the jet and the soft functions with 0 < a < 2, the coefficients are given a 
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At NLL order, the quantities in Eq. (88) are written as 
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and fi) is obtained from Ai{fii, fi) when aj are replaced by bj. 
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